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PASSIVE  INFRARED  SURVEILLANCE 
PART  I:  MODEL  FORMULATION 


INTRODUCTION 

This  paper  describes  an  analytical  model  for  predicting  the  performance  of  a  particular 
class  of  infrared  sensors  generically  described  as  infrared  search  and  track  (IRST)  devices.  An 
IRST  system  generally  consists  of  one  or  more  photodetectors  located  in  the  focal  plane  of  a 
scanning  optical  telescope,  and  a  complement  of  signal  processing  electionics  to  process  the 
detected  photocurrents.  The  signal  processor’s  task  is  to  determine  whether  or  r  ot  an  object  of 
a  particular  type  (a  "target")  is  anywhere  in  the  sensor's  field  of  view,  while  keeping  the  fre¬ 
quency  of  false  target  reports  to  an  acceptably  low  level. 

Thus,  the  Neyman-Pearson  criteria  conventionally  applied  to  radar  systems  is  appropriate 
also  as  an  objective  of  IRST  processor  design:  the  target  detection  probability  ( PD )  should  be 
maximized  for  a  giver,  maximum  tolerable  false  alarm  rate  (FAR). 

The  IRST  is  a  nonimaging  device,  as  contrasted  with  forward  looking  infrared  (FLIR) 
imaging  systems  {1,2].  The  search  and  track  device  may  be  required  to  keep  a  full  hemisphere 
(2ir  steradians)  under  constant  observation,  to  have  r  resolution  of  one  milliradian  or  less,  and 
to  operate  without  human  assistance  for  long  periods  of  time.  The  challenge  this  presents  to 
the  system  designer  is  further  magnified  by  the  abundant  opportunities  for  target/background 
confusion  offered  by  such  typical  background  scenes  as  cloudy  skies  and  cities. 

Previous  attempts  to  model  background  effects  on  IRST  system  performance  have  focused 
on  the  Wiener  spectrum  approach  (3-9],  a  frequency  domain  technique  originally  developed  for 
calculating  the  noise  variance  in  communication  circuits.  Unfortunately,  unlike  the  noise 
processes  typically  assumed  in  statistical  communication  theory  [10,11],  the  IRST  photocurrent 
is  a  highly  non-stationary  random  process.  As  the  sensor  scans  across  a  structured  background, 
the  spatially  non-uniform  scene  brightness  is  mapped  into  a  photocurrent  whose  mean  and  vari¬ 
ance  are  both  functions  of  time.  The  inadequacy  of  the  Wiener  spectrum  method  under  thetse 
conditions  has  been  appreciated  for  many  years  (12,13). 

The  IRST  model  developed  in  this  paper  requires  as  input  complete  descriptions  of  the 
IRST  sensor  and  the  scene  radiance  distribution  (possibly  including  a  target).  As  output,  the 
model  generates  the  probability  that  the  IRST  device  declares  (rightly  or  wrongly)  a  target’s 
presence  in  the  scene.  A  priori  knowledge  as  to  whether  a  target  was  in  fact  present  in  the 
specified  scene  allows  interpretation  of  the  probability  of  target  declaration  either  as  a  probabil¬ 
ity  of  target  detection  ( P„ )  or  as  a  "false  alarm." 

Photon  fluctuation  noise  is  the  only  stochastic  aspect  of  the  model;  the  background  must 
be  specified  as  a  radiance  map  of  arbitrary,  but  deterministic,  structure.  Thus,  the  performance 
predictions  made  with  this  model  are  background-conditional. 
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Previous  works  describing  analysis  techniques  for  IRST  systems  (e.g.,  Refs.  14  and  IS) 
are  restricted  in  validity  to  uniform  scenes.  The  contribution  of  this  paper  is  an  original 
method  for  calculating  IRST  performance  (i.e.,  the  parameters  PD  and  FAR)  that  is  inherently 
applicable  to  non-uniform  scenes.  The  method  described  here  can  be  used  to  assess  the  rela¬ 
tive  merits  of  a  variety  of  IRST  system  concepts,  all  operating  against  a  particular  infrared 
scene.  Alternatively,  an  IRST  sensor's  background-conditional  performance  can  be  evaluated 
for  each  member  of  an  ensemble  of  scenes,  in  order  to  establish  ensemble  average  performance 
statistics  for  a  realistic  range  of  operational  environments. 

Before  launching  into  the  main  body  of  analysis,  a  treatment  of  the  elementary  concepts 
involved  is  first  presented  in  the  next  section.  The  latter  part  of  the  subsection  (Current  Statis¬ 
tics)  presents  a  brief  discussion  of  how  values  for  the  mean  and  variance  of  a  photocurrent 
(i.e.,  the  "current  statistics")  are  obtained  from  knowledge  of  the  brightness  of  a  presumably 
uniform  background  scene. 

The  final  part  of  the  next  subsection  (Crossing  Rates  for  Fixed  Threshold  Detection) 
presents  Rice's  well-known  equation,  which  is  Eq.  (9)  in  this  text,  for  the  threshold  crossing 
rate  of  a  stationary  Gaussian  random  process  (16],  Evaluation  of  Eq.  (9)  requires  knowledge  of 
the  mean  value  and  variance  of  the  Gaussian  process.  A  description  of  IRST  performance 
against  uniform  background  scenes  is  obtained  by  inserting  the  current  statistics  from  the 
Current  Statistics  Subsection  into  Rice’s  equation,  Eq.  (9), 

Rice's  equation,  by  itself,  is  devoid  of  physical  content:  it  applies  equally  well  to  any  sta¬ 
tionary  Gaussian  process.  All  of  the  physical  parameters- the  optical  and  electrical  characteris¬ 
tics  of  the  IRST  sensor  and  the  radiance  of  the  background  scene— are  introduced  through  Eqs. 
(3)-(S)  for  the  current  statistics. 

The  method  described  above,  using  Rice's  equation  to  evaluate  search  set  performance 
against  uniform  scenes,  is  well  known  (4,14),  Since  a  target’s  presence  in  the  scene  would 
necessarily  render  the  scene  non-uniform,  and  since  the  method  based  on  Rice's  equation  is 
valid  only  for  uniform  scenes.  Rice's  equation  can  be  used  to  calculate  FAR,  but  not  PD. 

The  first  method  presented  for  calculating  search  set  performance  against  non-uniform 
scenes  is  gi\en  in  the  Non-Uniform  Scenes  Subsection.  Although  this  method  is  a  simple, 
heuristic,  extension  of  the  well-known  uniform  background  result  reviewed  in  the  subsection 
on  Uniform  Scenes,  it  appears  to  be  original.  The  heuristic  approach  of  the  subsection  on 
Non-Uniform  Scenes  has  the  advantage  of  being  both  easy  to  understand  and  easy  to  apply. 
Moreover,  comparison  with  rigorous  methods  indicates  that  the  heuristic  method  yields  numeri¬ 
cally  accurate  performance  predictions  as  long  as  the  variations  in  the  background  scene  are  not 
too  rapid  (17,18]. 

The  IRST  performance  model  developed  in  the  Analysis  Section  is  a  rigorous  generaliza¬ 
tion  of  the  Rice  equation  method.  In  fact,  a  point-by-point  correspondence  can  be  established 
between  the  simple,  well-known  analysis  of  the  Uniform  Scenes  Subsection  and  the  more  gen¬ 
eral,  new  analysis  in  this  report. 
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Appendix  A  generalizes  Eq.  (3). 

Appendix  B  generalizes  Eqs.  (4)  and  (S). 

Eq.  (55),  originally  derived  by  Cramdr  and  Leadbetter  (19),  generalizes  Eq.  (9).: 

As  expected,  the  analysis  in  this  report  redices  exactly  to  the  simple  Rice  equation 
method  in  the  uniform  background  limit.  More  generally,  the  method  developed  in  this  report 
is  also  applicable  to  non-uniform  scenes  (i.e.,  scenes  containing  targets  and/or  structured  back¬ 
grounds), 

The  techniques  developed  in  this  report  are  applied  to  the  evaluation  of  a  number  of  can¬ 
didate  signal  processing  structures  in  a  companion  paper  to  the  present  work  (18). 

It  is  noted  that  tracking  algorithms  are  neither  modeled  nor  discussed,  although  they  may 
play  an  important  role  in  clutter  rejection  as  well  as  target  tracking  (20),  Moreover,  no  con¬ 
sideration  is  given  to  the  availability  of  target  and  backgound  infrared  radiance  data  suitable  as 
model  inputs.  No  strategies  ire  proposed  for  synthesizing  1RST  processor  structures  to  satisfy 
either  a  priori  perform  once  requirements  or  optimization  criteria.  The  contribution  of  this  paper 
is  the  formulation  of  a  performance  analysis  model:  Complete  descriptions  of  an  IRST  sensor 
and  the  radiance  distribution  of  a  particular  scene  are  required  as  inputs.  As  output,  the  model 
generates  the  probability  that  the  IRST  device  declares  (rightly  or  wrongly)  a  target's  presence 
in  the  given  scene. 

SCANNING  BLIP  SENSORS;  ELEMENTARY  CONCEPTS 
A  Basle  Thresheld  Receiver 

In  order  to  provide  a  frame  of  reference  for  the  following  discussion,  it  is  necessary  to 
describe  a  simple  IRST  receiver  structure,  and  to  define  the  parameters  used  to  characterize 
IRST  performance. 

The  probability  that  the  IRST  device  makes  a  target  declaration  when  a  target  is  in  fact  in 
the  sensor's  Aeld  of  view  is  called  the  Probability  of  Detection  (P,>).  The  average  rate  at  which 
false  target  declarations  occur  is  termed  the  false-alarm  rate  (FAR). 

A  basic  threshold  comparison  receiver  is  shown  in  Fig.  1,  The  current  X(i)  at  the  output 
of  the  detector  is  input  to  an  electrical  Alter  of  transfer  function  H(f).  The  output  current 
Y(i)  of  the  electrical  Alter  is  compared  with  a  threshold  Y„(i).  If  Y(t)  exceeds  the  threshold, 
the  presence  of  a  "target''  is  declared;  otherwise,  no  target  declaration  is  made.  Target  detec¬ 
tions  and  false  alarms  are  both  manifested  as  threshold  crossings,  suggesting  the  following 
approach  to  IRST  performance  assessment: 

•  The  expected  number  of  threshold  crossings  that  a  particular  processor  experi¬ 
ences  against  a  given  infrared  scene  is  Arst  calculated.  As  discussed  in  Appendix 
C,  the  expected  number  of  threshold  crossings  during  a  given  time  interval  may 
be  interpreted  as  a  probability  of  target  declaration. 
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Fi|  I  -  A  basic  threshold  comparison  receiver.  The  photodetecior  in  this 
figure  is  ‘idealized’,  in  the  sense  that  it  is  presumed  to  have  a  perfect  all-pass 
electrical  frequency  characteristic;  the  frequency-dependent  part  of  the  detector 
responsivity  is  lumped  together  with  the  transfer  function  of  the  post-detector 
filter  to  obtain  HO  A  ‘target  declaration*  is  made  whenever  the  filtered 
current  Ki)  esceeds  the  threshold  level  K„(<) 

•  If  the  specified  scene  is  known  to  contain  a  target,  the  computed  probability  of 
target  declaration  is  interpreted  as  a  Probability  of  Detection  (/»/,);  otherwise,  a 
False  Alarm  Rate  (FAR)  interpretation  is  given. 

As  will  presently  be  discussed,  it  is  highly  desirable  that  the  threshold-establishing 
mechanism  suppress  clutter-induced  threshold  crossings  by  increasing  Y„(i )  when  Y(t)  is 
"clutter-like." 

Rather  than  allow  Y„(i)  to  take  on  an  a  priori  constant  or  functional  value,  it  is  necessary 
to  establish  the  threshold  by  some  means  that  "adapts"  K„(r)  to  the  prevailing  background  con¬ 
ditions. 

A  similar  type  c<f  signal  processing  problem  has  been  addressed  in  the  radar  [21,22]  and 
sonar  [23]  literatures.  A  candidate  adaptive  threshold  scheme  adapted  from  the  earlier  work 
[24]  is  deputed  in  Fig.  2. 

The  block  with  transfer  function  exp  (- j2wfTd )  introduces  a  delay  of  Td  seconds.  The 
triangular-shaped  block  in  this  figure  denotes  an  ideal  all-pass  amplifier  of  gain  K.  Note  that 
the  delay  time  Td ,  gain  A,  and  transfer  functions  //„(/)  and  //(/)  are  all  design  variables. 
Strategies  for  choosing  the  design  variables  in  order  to  satisfy  particular  performance  require¬ 
ments  or  optimization  criteria  will  not  be  discussed.  This  report  is  devoted  to  developing  a  for¬ 
mulation  for  the  expected  number  of  threshold  crossings  for  IRST  receivers  structured  as  in 
Fig.  2,  under  the  assumption  that  the  design  variables  have  .ill  been  specified. 
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Fig  2  —  A  simple  adaptive  threshold  scheme  is  illustrated  The  transfer 
function  espl  /2ir/  7]/)  introduces  a  delay  or  t],  secant  ensuring 
decorrelation  of  the  random  processes  Hi)  and  1 1 )  The*  uticunce  of 
the  random  currents  is  seen  by  inspection  of  Fig  I 


Uniform  Scenes 


Current  Statistics 

Whe”  the  sensor  of  Fig.  1  scans  across  a  uniform  background  scene,  the  output  current 
Y(t)  is  “statistically  stationary.”  The  meaning  of  statistical  slalionarity  will  now  be  discussed  us 
background  to  the  discussion  of  non-uniform  scenes  and  non-stationary  processes  thut  follow  in 
the  Non-Uniform  Scenes  subsection.; 

It  is  assumed  thut  the  sensor  is  scanned  arid  re-scanned  over  the  same  scene,  and  that 
there  arc  no  chunges  in  cither  the  scene  or  the  sensor  from  one  scan  to  the  next. 

The  current  .»•(/)  during  the  course  of  any  one  particular  scan  is  culled  a  "sample  func¬ 
tion"*  of  the  random  process  Yti)  (The  process  )'(/),  in  turn,  may  be  thought  of  us  the 
infinite  ensemble  of  possible  sample  functions.)  The  current  sumplc  function  obtained  on  the 
n"'  scan  is  designated  .►’(/;«).  We  now  consider  a  particular  one  of  these  sumplc  functions, 
»  </  P,  depicted  in  Fig.  3. 


The  time  variations  in  i  (/;l)  have  their  origin  in  the  time-of-urrivul  fluctuations  of  the 
individual  photons  incident  on  the  detector. t  Thus,  the  fluctuations  in  t(/,P  ure  independent 
of  the  scan  velocity  and  are  present  regardless  of  whether  the  sensor  is  scanning  or  motionless 
The  average  current  at  a  particular  instant  of  lime  /„  may  be  defined  us  the  "ensemble  average": 


/«,(/„)=  tim  -!  =flf(/„)l 

N  N  ,7?\ 


In  order  for  Y(t)  to  be  a  stationary  process,  A  is  necessary  that  m,tr),  as  defined  in  Kq.  (I),  be 
independent  of  time.  Thus, 

in, tl„)  -  /«,(/().; 


Consistent  with  j  vonvenliott  ol  random  process  theorv ,  sIihIusik  quantities  an-  assigned  lapilul  lelteis,  m  it  It  sample 
values  designated  bv  the  corresponding  lower  vase  letters 

*|his  kind  ol  noise  is  oltet'  called  'photon  lliiituation  noise"  or  "quantum  noise"  When  the  quaiium  noise  assiKialed 
with  the  background  light  is  the  dominant  n.nse  tvpe  in  the  seiism..  the  sensor  is  said  to  be  operating  in  the  "B.uk 
ground  l  muted  Perlorniancc*  tBI  lt’>  regime 
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Fi|.  3  -  Three  umple  functions  of  the  random 
current  Y(i)  are  shown  as  functions  of  time.  The 
current  Y(i)  is  taken  at  the  output  of  the  post¬ 
detector  Alter  as  the  sensor  is  scanned  over  a  uniform 
scene  (cf.  Fig.  I).  The  sample  functions  are  desig¬ 
nated  v(r;»).  «- 1,2.3.  .  ..  The  sample  functions 
display  random  time  variations  caused  by  time-of- 
arrival  fluctuations  of  the  individual  photons  incident 
on  the  detector.  Since  the  current  YU)  is  stationary, 
the  ensemble  average  mean  value  (variance)  defined 
by  Eq.  (I)  (Eq.  (2))  is  the  same  at  lime  /„  as  at  time 
where  times  i„  and  i\  are  arbitrary. 


where  the  times  i„  and  /,  are  totally  arbitrary  (cf.  Fig.  2),  Similarly,  the  mean-square  deviation 
of  Y(t)  from  its  average  value  (i.e„  the  "variance"  of  2)  may  be  defined  at  each  instant  of  time 
as 


«■{(/)  =  lim 
,v-<» 


1 

N 


,Y 


I 


l  YU\n )  -  my(/»* 


=  £{(K(f)-/My(/)|J). 


(2) 


The  variance  <r?,  like  the  mean  m>,  is  independent  of  time  for  stationary  processes. 

Equations  (I)  and  (2)  are  satisfactory  for  illustrating  the  concept  of  "ensemble  averaging"; 
however,  it  is  desireable  to  have  a  different  means  for  actually  calculating  the  values  of  w>  and 
<t  j  in  terms  of  standard  background  and  sensor  parameters. 

It  is  first  necessary  to  define  the  average  value  my  of  the  current  X(i)  (cf.  Fig.  I): 

my  -  £{<¥(/)}  -  rje/wy,  (3) 


where 


V 

e 


—  quantum  efficiency  of  the  detector 


—  electronic  charge 


coulombs 

electron 


electrons 

photon 
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and 

m0-  average  background  photon  flux  ^.?ion.s.  incident  on  the  detector. 

see 

It  may  be  shown  that  for  sensors  operating  in  the  BLIP  regime  the  mean  value  and  variance  of 
Yit),  originally  defined  by  Eqs.  (1)  and  (2),  can  be  calculated  in  terms  of  mx  as  follows: 

mY-HiQ)mx.  (4) 

and  (251 

<r\-lemxtif,  (5) 

where  //( 0)  is  the  zero-ordinate  of  the  transfer  function  Hif)  (cf.  Fig.  1),  and  A/ is  the  noise 
bandwidth  of  //(/).  For  bandpass  Hi/),  Hi 0)  -  0.  It  follows  from  Eq.  (4)  that  mY  -  0  for 
this  case. 

Since  the  scene  is  spatially  uniform,  the  average  photon  flux  mQ  is  independent  of  time. 
It  follows  from  Eqs.  (3)  •  (3)  that  the  mean  mY  and  variance  are  also  independent  of  time, 
justifying  the  claim  of  stationarity  for  the  current  YU). 

Assuming  that  the  transfer  function  H  if)  is  normalized  as  follows: 

max  Hif)  -  1.  (6) 

the  noise  bandwidth  A/ in  Eq.  (S)  may  be  calculated  from  the  equation: 

(7) 

Crossing  Rales  for  Fixed  Threshold  Detection 

It  is  now  assumed  that  the  signal  processor  of  Fig.  1  is  implemented  such  that  the  thres¬ 
hold  y„  is  equal  to  a  constant.  The  fixed  threshold  y„  is  depicted  on  the  sample  function  plots 
of  Fig.  3.  A  brief  outline  will  now  be  given  of  a  method  for  calculating  the  average  number  of 
times  ntj  that  the  random  process  YU)  crosses  the  threshold  during  a  time  interval  of  duration 
T  seconds.  (The  relationship  between  the  mean  number  of  crossings  mj  and  the  usual  search 
set  performance  parameters  PD  and  FAR  is  discussed  in  Appendix  C.) 

The  expected  number  of  threshold  crossings  mjiO.T)  during  the  time  interval  |r|<r/2 
may  be  written  in  terms  of  a  'crossing  rate"  mj  as: 

mjiO.T)  -  thjT  (8) 

According  to  Rice  (16),  thj  may  be  calculated  as 

t”j  “  i»jo  exp  i-u1/ 2),  (9) 

where 
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with  A/ given  by  Eq.  (7).  Also, 

u  =  (y0~nty)/<r  y,  (11) 

with  ttt  y  and  <r  y  given  by  Eqi>.  (4)  and  (S). 

The  quantity  u  defined  by  Eq.  (11)  may  be  thought  of  as  a  normalized  threshold  level. 

It  follows  from  Eq.  (9)  that  the  expected  number  of  threshold  crossings  drops  off  rapidly 
as  the  theshold  level  is  increased 

Finally,  it  should  be  noted  that  the  uniform-background  threshold  crossing  formalism 
described  above  cannot  be  used  to  calculate  the  expected  number  of  threshold  crossings  associ¬ 
ated  with  a  target’s  presence  in  the  scene*,  because  a  target's  presence  would  render  the  scene 
non-uniform. 

Non-llniform  Scenes 

Current  Statistics 

It  is  now  assumed  that  the  sensor  of  Fig.  1  is  scanned  a  number  of  times  over  the  same 
non-untform  scene,  and  that  there  are  no  changes  in  either  the  scene  or  the  sensor  from  one 
scan  to  the  next.  A  number  of  sample  functions  of  the  resulting  current  process  YU)  are  dep¬ 
icted  in  Fig.  4. 

Once  again,  the  ensemble  average  mean  and  variance  of  Y(t)  are  defined  by  Eqs.  (1)  and 
(2),  However,  as  discussed  next,  Y(t)  is  now  a  non-stationary  process,  i.e.,  my  and  are 
functions  of  time. 

As  an  illustration  of  how  such  non-stationary  processes  arise,  it  is  now  assumed  that  the 
infrared  scene  encompasses  regions  of  blue  rky  and  clouds  where 

">0('av)"*  average  photon  flux  incident  on  the  detector  when  blue  sky  is  being  observed, 
and 

myO,)-  average  photon  flux  incident  on  the  detector  when  cloud  is  being  observed. 

The  photon  flux  wy  in  Eq.  (3)  is  seen  to  be  a  function  of  time:  /wy(/)  takes  on  the  value 
mQUl )  at  a  time  /(  when  a  cloud  is  in  the  field  of  view,  and  it  takes  on  a  different  value 
myO/jv)  at  a  time  f#s  when  the  scanning  Add  of  view  includes  only  blue  sky.  Thus,  the  process 
YU)  is  non-stationary  when  the  scene  is  non-uniform,  because  the  mean  and  variance  of  Y(t) 
are  seen  from  Eqs.  (3)-(5)  to  be  functions  of  timet. 

As  discussed  in  Appendix  C,  the  incremental  number  of  threshold  crossings  associated  with  a  target’s  presence  in  the 
scene  provides  an  estimate  of  the  conventional  search  set  parameter  PD  By  definition,  P0  is  the  probability  of  larget 
detection 

Vor  the  present,  it  suffices  lo  say  that  the  forms  of  Eqs  (41  and  (5)  indicate  that  a  time-varying  must  give  rise  lo 
time  varying  mr  and  <ry  However,  it  should  be  noted  that  Eqs  (4)  and  (S)  arc  only  strictly  valid  for  stationary 
processes,  le,  for  time-invariant  mx  Generalizations  of  Eqs  (4)  and  (5)  valid  for  both  stationary  and  non-siaiionary 
processes  are  given  by  Eqs  (B-30)  and  (B-31) 
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Fig.  4  —  Three  simple  functions  of  the  ran¬ 
dom  current  Y(i)  are  shown  as  functions  of 
time.  This  figure  is  similar  to  Fig.  3.  except 
the  sensor  is  now  presumably  scanned  over 
a  non-uniform  scene  The  time-varying 
ensemble  average  mj.fr)  of  Yfi)  is  shown  as 
a  dished  curve  superposed  on  each  of  the 
three  depicted  sample  functions  (solid 
curves).  The  ensemble  average  my  is  still 
defined  by  Eq.  (I);  however,  the  fact  that 
my  is  now  a  function  of  time  implies  that 
)'(()  is  now  a  nonstationary  random  pro¬ 
cess. 


The  time-varying  mean  value  nty(t)  is  superposed  as  a  dashed  curve  on  each  of  the  sam¬ 
ple  functions  >>(/)  depicted  in  Fig.  4. 

Crossing  Kales  for  Fixed  Threshold  Detection 

The  performance  of  a  fixed  threshold  signal  processor  (cf.  Fig,  1)  against  a  non-uniform 
scene  can  be  characterized  in  terms  of  the  quantity  ntj ,  where* 

mj(O.T)  =£UI=  /,  mj(t)dr.  (12) 

where  T„  is  the  time  interval  |/|  <  772,  J  is  an  integer  random  variable  equal  to  the  number  of 
times  that  the  current  Y(t)  crosses  the  threshold  level  ,v„  during  the  time  interval  Tin  and  £’(.) 
is  the  statistical  expectation  operator  as  defined  in  Eqs.  (I)  and  (2).  Equation  f  1 2)  is  a  straight¬ 
forward  generalization  of  Eq.  (8)  to  allow  for  the  possibility  of  time-variable  threshold  crossing 
rates  thj. 

As  long  as  the  time  variation  of  nix(t)  is  slow  compared  to  the  time  variation  of  the 
impulse  response  li(t)  of  the  post-detector  filter  (cf.  Fig.  1),  a  good  estimate  for  hij(t)  can  be 
obtained  from  Eq.  (9). 

Cl  Appendix  C  for  j  discussion  ol  the  relationship  ul  nijU>.  1)  to  the  usual  IKS  f  performance  parameters  I’/,  and 

FAR 
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The  following  steps  are  then  followed  in  calculating  mj(0,T);. 

a.  The  time-varying  mean  current  mx{t)  is  derived  from  the  time-varying  photon  irradi- 
ance  mQ(i)  by  means  of  Eq.  (3).:  (A  detailed  formulation  for  mxU)  in  terms  of  the  back¬ 
ground  radiance  distribution  is  provided  in  Appendix  A.) 

b.  Estimates  of  mr(t)  and  <r7y(t)  are  obtained  from  Eqs.  (4)  and  (5).:  (More  rigorously, 
Eqs.  (B-30)  and  (B-31)  may  be  used  to  obtain  my(t)  and  o-f(r).) 

c.  Equation  (11)  is  evaluated  for  the  time-varying  normalized  threshold  u(t).. 

d.  Equation  (9)  is  evaluated  for  the  time-varying  threshold  crossing  rate  rrij(t) , 

e.  Equation  (12)  is  evaluated  for  the  expected  number  of  threshold  crossings  w7(0,  T)  , 

Numerical  examples  following  the  above  prescription  typically  show  that  the  crossing  rate 
function  rhj(t )  is  extremely  sharply  peaked  (cf.  Fig.  S),  Consequently,  appreciable  contribu¬ 
tions  to  mjiO.T)  only  accrue  in  the  near  neighborhood  of  points  such  as  tp  in  Fig.  S.  It  is 
shown  in  Ref,  17  that  the  time  tp  in  Fig.  Sb  is  a  saddle  point  of  the  crossing  rate  integral  Eq. 
(12),  and  that  Eq.  (12)  may  be  approximated  asymptotically  as: 

mj( 0.  T)  =*  ntjitp)  8/p,  (13) 

with  rrijii,,)  obtained  from  Eq.  (9).;  The  quantity  btp  is  the  effective  interval  of  time  during 
which  mY(i)  remains  in  the  near  neighborhood  of  its  peak  value,  from  the  standpoint  of  cross¬ 
ing  rate  calculations.  An  expression  for  htp  is  derived  in  Ref.:  17. 

The  implications  of  Eq.  (13)  for  system  performance  are  illustrated  with  the  aid  of  Fig.  6. 


Fig  5  —  Pari  a)  is  an  illusiralive  plot  ol'  niy(t)  vs  /, 
where  is  the  mean  value  of  the  tillered  current 
>'(/)  Also  shown  is  a  constant  threshold  current 
lying  above  the  peak  value  of  my  The  function 
myil)  takes  on  Us  peak  value  at  the  time  ip 

Part  (b)  is  a  plot  ol  the  threshold  crossing  rate 
in  jit)  corresponding  to  the  threshold  v„  and  mean 
current  wy(i)  of  part  (a).  The  entire  contribution 
to  the  crossing  rale  integral,  F.q  (12),  accrues  in  tlie 
very  near  neighborhood  of  ip 
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Pig.  6  —  This  figure  illusirates  a  critical  shortcoming  of  constant 
threshold  processing.  The  slowly  varying  maximum  centered  at  in 
presumably  has  its  origin  in  the  non-uniform  background  scene.  The 
narrower,  lower  amplitude  spike  centered  at  i„  is  due  to  a  target. 

The  likelihood  of  a  false  alarm  (i  e..  a  clutter-induced  threshold  cross¬ 
ing)  grows  rapidly  as  the  threshold  level  v„  is  reduced  There  is  no 
way  for  the  constant  threshold  processor  to  detect  the  target  peak  at 
i„  without  also  incurring  a  false  alarm  arising  from  the  clutter  peak 
centered  at  r ,, 

A  plot  of  the  threshold  crossing  rate  iitj(i)  corresponding  to 
this  figure  would  show  that  the  probability  of  a  threshold  crossing, 
and  hence  a  false  alarm,  is  far  greater  at  time  than  at  any  other 
time 

The  large,  relatively  slowly  varying  maximum  centered  at  in  Fig.  6  is  presumed  to  have 
its  origin  in  the  background  scene.  The  narrower,  lower  amplitude  spike  centered  at  tn  in  Fig.  6 
is  presumed  to  be  due  to  a  "target. " 

It  follows  from  Eqs.  (13)  and  (9)  that  the  likelihood  of  a  clutter-induced  threshold  cross¬ 
ing  grows  rapidly  as  the  threshold  level  y„  in  Fig,  6  is  lowered.  A  clutter-induced  threshold 
crossing  (i.e. ,  a  "false  alarm")  becomes  a  virtual  certainty6  when  there  is  a  "mean-crossing,"  i.e., 
when  the  threshold  level  actually  intercepts  the  mean  current  m,(t).  117,18)  There  is 
apparently  no  way  for  the  constant  threshold  processor  to  detect  the  target  peak  at  t„  without 
also  incurring  a  false  alarm  arising  from  the  clutter  peak  centered  at  in , 

Crossing  Rales  for  Constant  False  Alarm  Rale  (CFAR)  Adaptive  Threshold  Detection 

The  performance  of  an  adaptive-threshold  processor  is  illustrated  with  the  aid  of  Fig.  7. 

The  processor  is  presumed  to  have  some  means  for  deriving  high-confidence  estimates  for 
/»,(/)  and  <r  > ( r ) ,  defined  as  and  <r>,  respectively.  When  m>(/)  is  "slowly-varying”  the  pro¬ 
cessor  establishes  yjt)  as: 


y„(t)  =  /«>(/)  +  Ktr ,(/). 


(14) 


* 

This  has  been  established  by  integrating  Ihe  crossing  rate  I  unction  nij  over  an  interval  ot  time  containing  a  lime  point 
im  tor  which  i ii)(im  )  -  i„  An  asymptotic  analysis  (cl  Ref  17)  has  shown  that  he  crossing  rale  integral  I  q  (12)  is 
incremented  by  unity  for  each  such  lime  i„„  contained  in  the  interval  ol  integration  /„  I  his  result  is  not  at  all  surpris¬ 
ing,  and  may  '»e  taken  as  evidence  that  the  theory  developed  here  is  consistent  with  common  sense 
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my(t) 


*0 

Fig.  7  -  This  figure  illustrates  an  important  potential 
advantage  of  adaptive  threshold  processing.  The  mean 
current  m^l/l  is  the  same  as  for  Fig.  6,  with  a  clutter  peak 
centered  at  time  and  a  lower-amplitude  target  spike  cen¬ 
tered  at  time  The  adaptive  threshold  i„(/)  is  presum¬ 
ably  able  to  accurately  track  the  slowly  varying  background 
signal,  but  not  the  more  rapidly  varying  target  signal. 

Thus,  target  detection  is  assured,  while  the  probability  of  u 
false  ularm  is  kept  acceptably  small  As  contrasted  with  the 
situation  of  Fig.  6,  the  probability  of  a  false  alarm  is  now 
no  greater  in  the  neighborhood  of  time  than  at  any  other 
lime. 

The  threshold  is  interpolated  through  periods  of  "rapidly  varying"  rn>(/)  by  means  of  a  smooth¬ 
ing  filter  (cf.,  for  example,  the  neighborhood  of  /„  in  Fig.  7).  The  adaptive-threshold  constant 
K  in  Eq.  (14)  is  a  design  parameter. 

From  Eqs.  (11)  and  (14), 

«(/)  =»  K[(r  >0)/»r  >(f)|  +  |>«,0)  -  /Wj(f)j/rr  ,(f).  (15) 

When  the  estimation  errors  are  sufficiently  small. 


and 


(>,(/)  =  (TyU).  (16) 

/«>(/)  =S  /« y  ( r ) . 


it  follows  from  Eq.  (15)  that 


uU)  —  K. 


(17) 


Thus,  «(/)  is  rendered  time-invariant  by  the  adaptive-threshold  processor  when  there  is  no 
rapidly- varying  target  contribution  to  /H»(r),  When  there  are  no  targets  in  the  scene,  the  mean 
current  /w » ( / )  is  assumed  to  be  slowly-varying,  and  the  expected  number  of  crossings  during 
the  time  interval  T„  may  be  calculated  from  Eqs.  (9),  (12),  and  (17)  as 

ntjiO,  T)  =  J  ritj„  exp  (-K1/ 2)  ill.  (18) 

Since  both  mJo  and  K  are  time-invariant,  Eq.  (18)  becomes 

'"A 0.  T)  =  ihj„  exp  (-K2/ 2)  T.  (19) 

Equation  (19)  has  the  same  form  as  the  crossing  rate  expression  for  uniform  scenes,  Eq.  (8). 
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The  kind  of  processor  just  described  has  been  called  a  constant  false  alarm  rate  (CFAR) 
processor,  since  the  threshold  crossing  rate  is  now  independent  of  time,  i.e.,  a  crossing  is  no 
more  likely  to  occur  when  scanning  a  region  of  non-uniform  background  than  when  scanning  a 
region  of  uniform  background.  For  example,  with  reference  to  Fig.  7,  the  crossing  rate  is  now 
no  greater  at  ip  than  at  any  other  time. 

Unfortunately,  the  CFAR  processor  is  generally  a  non-realizable  ideal:  it  has  been 
assumed  that  the  processor  is  able  to  estimate  the  quantities  mYU)  and  <rYU)  to  as  high  a  pre¬ 
cision  as  desired.  Errors  in  the  estimated  values  for  mY  and  <r  Y  are  usually  unavoidable,  giving 
rise  to  appreciable  time-dependence  in  Eq.  (IS)  for  uU). 

ANALYSIS 

Introduction 

The  objective  of  the  next  subsection  is  to  present  an  expression  for  the  expected  value  of 
the  number  of  times  a  nonstationary  noise  YU)  crosses  a  nonstationary  stochastic  threshold 

W'), 

The  crossing-rate  formulation  that  results,  Eqs.  (25)  and  (27),  requires  knowledge  of  the 
joint  density  function  fyyy  y  of  the  current  YU),  its  time  derivative  YU),  the  threshold  Y„U), 

\ 1  1 »  '  II 

and  its  time  derivative  Y0U). 

As  discussed  in  Ref.  17.  the  currents  YU)  and  Y0U)  are  non-stationary  Gaussian 
processes.  It  follows  that  YU),  YU),  Y0(r),  and  Y0U),  are  jointly  Gaussian  processes.  The 
joint  density  fyyy .  ^  can  thus  be  expressed  in  terms  of  a  covariance  matrix  A  .< 

Assuming  that  the  Altered  current  YU)  and  the  threshold  Y0U)  are  uncorrelated 
processes,  the  fourth-order  density  function  factorizes  into 

f Wr„  ir„  “  fry  fy0  r/  ^0) 

The  justiAcation  for  Eq.  (20)  is  discussed  in  the  next  section.  Expression  (S2)  for  the  jointly 
Gaussian  fYY  is  then  used  with  the  general  crossing  rate  Expression  (48)  to  derive  a  more 
explicit  crossing  rate  expression,  Eq.  (55).,  Evaluation  of  Eq.  (55)  for  mj  requires  the  expres¬ 
sions  derived  in  Appendix  B  for  the  time-varying  current  statistics  <rzU).  <jyU),  and  rU)  (cf. 
Eqs.  (B-30)  -  (B-34)).; 

The  complete  expression  for  the  average  crossing  rate,  Eq.  (47),  generally  requires  the 
numerical  integration  of  a  somewhat  complicated  integrand.  The  section  called  Crossing  Rates 
for  Adaptive  Threshold  Processors  is  devoted  to  deriving  an  approximation  to  Eq.  (47),  The 
result,  Eq.  (67),  is  the  principal  analytical  result  of  this  report. 

A  Basic  Equation  for  Curve  Crossing  Rates 

The  integer  random  variable  J  is  deAned  as  the  number  of  zero-crossings  of  a  random 
process  GU)  on  a  time  interval  |r|  <  T/2.  The  expected  value  of  J  is  deAned  as  mj( 0,  T). 
Thus, 


mj( 0,  T)  -  E[JX, 
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where  £{.)  denotes  an  ensemble  average.  Defining  the  crossing  rate  function  mjit)  as  in  Eq. 
(12),  it  can  be  shown  that  (cf,  Ref.  26,  p.  514) 


-  J\g\fGa  (0 ,g)  dg , 


The  function  fGG  in  Eq.  (22)  is  the  joint  probability  density  function  of  the  process  Git)  and 
its  time  derivative  Git). 

Equation  (22)  is  well  known;  however,  most  references  to  it  appear  to  impose  a  stationer- 
ity  requirement  on  G  that  is  not  actually  necessary.  The  applicability  of  this  equation  to  nonsta¬ 
tionary  processes  appears  to  have  first  been  recognized  by  Cramdr  and  Leadbetter  [19]. 

The  domain  of  integration  in  Eq.  (22)  is  a  matter  of  some  interest.  If  one  wishes  to  cal¬ 
culate  only  the  expected  number  of  positive  slope  zero-crossings,  i.e.,  the  expected  number  of 
times  that  both 

Git,,)  -  0  (23) 


Git„)  >  0  (24) 

are  both  satisfied  on  the  interval  M  <  77 2,  the  lower  and  upper  limits  of  integration  in  Eq. 
(22)  should  be  chosen  as  0  and  «>,  respectively..  The  resulting  expression  for  mji 0,  T), 


mjiO.T)  -  Jr  dt  dg  UI/CG(0i)  . 


does  not  include  zero-crossings  of  the  type  depicted  in  Fig.  8,  for  which  g  <  0.  Apparently,  mj 
is  sensitive  only  to  the  "right  type”  of  zero  crossing,  as  defined  by  the  limits  of  integration  in 
Eq.  (22), 


Fig  8  —  A  sample  function  #(/)  of  the  random  process  GO)  is  dep¬ 
icted  as  a  function  of  lime.  The  particular  sample  function  chosen 
has  a  down-crossing  at  time  i„ 
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The  process  Git)  is  now  assumed  to  be  formed  as  the  difference  of  two  stochastic 
processes  Yit)  and  Y0it).  Thus, 

Git)  s  Yit)  -  Y0(t),  (26) 


where  Y0(t)  is  referred  to  as  the  threshold  process”.:  Without  making  any  assumptions  with 
respect  to  the  statistics  of  Yit)  and  Y0it)  (e.g.,  each  process  may  be  both  nonstationary  and 
non-Gaussian)  it  follows  from  Eq.  (26)  that 

/cdW-*)  ■  ffWl  frnarJC'S  +  n.(.v)-  (27) 


The  proof  of  this  equation  is  straight-forward  (cf.  Ref.  27,  p.  131).  Assuming  that  the 
processes,  Y.  Y,  Y0,  and  Y0,  are  jointly  Gaussian  (cf.  Ref.  17  for  justification),  their  joint  density 
can  be  expressed  in  terms  of  their  covariance  matrix  (27)  A.  The  matrix  A  has  four  rows  and 
four  columns,  for  a  total  of  sixteen  elements.  Written  in  partitioned  form,  we  have 


C  Cir 

c,  c 


where  the  superscript  T  denotes  the  matrix  transpose  operation.  The  submatrices  C,  C„,  and  Ct 
are  defined  as 


and 


C^(/,/)l 

Cyy(f»f  ) 

<r$0)  1' 

rl(t) 

Gy  y  it, I ) 

f  ii  *i# 

C  y  v  (/,/) 

'a 's 

<r\  it) 

C, 


Cf,  f(^)  Cyyit,t) 

Lyyit.t)  Gy  yit, t) 


(29) 

(30) 


(31) 


The  scalar  covariances  that  comprise  these  elements  are  defined  by 

C\p(/| ./*)  =  f|j^(/i)  -  «,«(/|)|  JflOj)  —  #w«(^z) j|. 

where  A  and  B  take  on  the  values  Y.  Y,  Yu,  and  Y„,  as  appropriate.  Also, 

mAit)  =  E{Ait)) 
and 

<rAit)  S  G^it.l).. 


(32) 

(33) 

(34) 


It  follows  from  Eq.  (32)  that  Cyyit.t)  -  CYyit,t)  and  Cr  Y  it.t)  m  CY  Y  it,t).  Thus,  the 
matrices  C  and  C„  are  symmetric. 
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The  elements  of  C  are  obtained  directly  from  Eqs.  (B-17),  (B-20),  and  (B-23): 

o-f(r)  -  e[mxU)  ®  A2(f)},  (35) 

<r \(t,t)  -  e[mxU)  ®  (Mr))2},  (36) 


and 


Cyy(t.l)  «  <Ty(t)  &YU), 


(37) 


where  e  is  the  electronic  charge,  hU)  is  the  impulse  response  of  the  post-detector  filter  (cf. 
Fig.  1),  mx(t)  is  the  mean  value  of  the  current  XU)  (cf.  Fig.  1),  and  ®  is  the  convolution 
operator: 

(37a) 


fU)  ®  git)  m  f  mfU-x)g(x)dx. 


The  relationship  of  mxU)  to  the  radiance  of  the  scene  under  observation  and  the  optical 
parameters  of  the  IRST  sensor  is  discussed  in  Appendix  A,  and  expressed  quantitatively  by  Eq. 
(A-2). 

Expressions  for  the  elements  of  Eq.  (30)  may  be  obtained  as  direct  adaptations  of  Eqs.  (35) 
(37).  It  follows  from  Fig.  2  that 

Y„U)  -  Kh„U~Tj)  ®  *(/),  (38) 

where  h„U)  is  the  Fourier  inverse  of  //„(/),  Taking  the  expected  value  of  both  sides  of  this 
equation,  we  have 

mYU)  m  Kh„U-Tj)  ®  mxU),  (39) 

which  is  analagous  to  Eq.  (BOO).  It  may  also  be  shown,  analagous  to  Eqs.  (35)  -  (37)  that: 

erf  U)  -  eKmxU)  ®  h}U-7d).  (40) 

<r f  U)  -  eKmxU)  OlKU-T,,))2,  (41) 

o 

Cr  Y  U.t)  m  <r  Y  U)  &Y  U).  (42) 


and 


It  remains  only  to  formulate  similar  expressions  for  the  elements  of  C(  in  order  to  complete  the 

specification  of  the  joint  density  fyXY  Y , 

*  *  '«  '« 

Assuming  that  the  processes  YU)  and  Y0U)  are  both  derived  from  the  process  XU)  by 
means  of  the  structure  shown  in  Fig.  2,  it  is  shown  in  Ref.  17  that,  if 


then 


KU-TJhU)  -0. 


C,  -0. 


(43) 


(44) 


That  is,  choosing  a  sufficiently  long  time  delay  Tj  in  Fig.  2  validates  the  factorization  of  the 
fourth  order  density  frrr  Y  into  the  product  of  two  second  order  densities: 


fyyy  y  CV>Ayo jo )  “  fYY  (yj)fY  y  U^ojo)’ 


(45) 
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From  Eqs.  (22),  (27),  and  (45), 


mj(t)  -  fYJn  ({,*))  ({,i+t))  . 


With  the  change  of  variable  z  -  g+t),  Eq.  (46)  may  be  written 


mj(t)  -  Ey  y  \mj(l k.j)0)|, 


where,  by  definition, 


J,oo 

„  \z-yn\fYY(yn.z)dZ. 

w» 

£^(.1  =  ffdidl)  fy  yJX.I))  {.) 


(46) 

(47) 

(48) 

(49) 


The  quantity  mj,  previously  interpreted  as  the  zero-crossing  rate  of  the  process  G  (cf.  Eq. 
(22)),  is  now  interpreted  as  the  threshold  crossing  rate  of  the  process  Y(t)  (cf.  Eq.  (47)). 

The  quantity 

mj(i\yu,y„)  =  mjU  I)  (50) 

defined  by  Eq.  (48)  will  be  referred  to  as  the  "threshold-conditional  crossing  rate." 
Non-Stationary  Gaussian  Processes 

Further  development  of  Eq.  (48)  for  the  threshold-conditional  crossing  rate  is  contingent 
on  obtaining  a  suitable  expression  for  the  joint  density  fYY.  The  objective  of  this  section  is  to 
evaluate  Eq.  (48)  for  the  particular  case  of  a  bivariate  Gaussian  density,  Eq.  (52). 

The  justification  for  assuming  a  Gaussian  distribution  for  jy  ('  nd  hence  for  fY  and  fYY  as 
well)  is  discussed  in  Ref.  17,  As  shown  in  Ref.;  17,  the  relative  error  in  ritj  is  approximately 
equal  to  the  relative  error  <  in  the  density  function  of  }'(/), 

fY{y)  2^-  (1  +  €).  (51) 

XT  , 

where  <f>(.)  is  the  Gaussian  density  function,  Eq.  (57),  The  Edgeworth  series  expansion  (16)  of 
fy(y)  provides  a  simple  and  easily  evaluated  expression  for  the  relative  error  c.  Sample  calcu¬ 
lations  described  in  Ref.  (17)  show  that  e  is  negligibly  small  for  typical  system  and  background 
parameter  values.  Thus,  the  joint  density  fYY  is  now  assumed  to  have  the  following  bivariate 
normal  form: 

/n(y».r)  *  2it(t  Y(t /(l— r2)l/J  exp  -  - - .  (52) 

2(1  -  r) 

where  the  quantities  mY,  <rK,  <rz ,  and  r,  are  obtained  from  Eqs.  (BOO)  -  (B04),  Also, 
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and 


U  =  (y,-my)/(T  y 


v  = 


Substituting  Eq.  (52)  into  Eq.  (48),  it  may  be  shown  that  (171 


(l-A2)172  4(u)  { <f>(p)  +  p  4> (/»)), 


(53) 

(54) 


(55) 


where,  by  definition: 

p  =  j(!-rJ)l/2  a, j  ( m ,  +  ru<rx-y„).  (56) 

The  functions  *(.)  and  <t»(.)  in  Eq.  (55)  are  defined  as  follows: 

4(x)  =  (2ir)~l/J  exp(— jt2/2)  (57) 

and 

*(x)  =  (58) 


The  crossing  rate  Eq.  (55)  is  originally  due  to  Cramlr  and  Leadbetter  (19). 

Equations  (47)  and  (55),  together,  represent  a  formal  means  for  calculating  the  mean  threshold 
crossing  performance  cf  the  adaptive  threshold  processor  depicted  in  Fig.  2.  However,  the 
evaluation  of  Eq.  (47)  appears  to  present  some  significant  calculations!  difficulties.  These 
difficulties  are  obviated  by  means  of  the  approximate  method  of  evaluation  pursued  in  the  fol¬ 
lowing  section. 

Crossing  Rates  for  Adaptive  Threshold  Processors 

Numerical  results  obtained  thus  far  indicate  that  the  correlation  coefficients  r  and  r„  are 
typically  much  less  than  unity,  as  follows: 

\r(l)\  -  |Cn(/./)  |<r  y(l)cr  y(/)l~'|  «  I  (59) 

and 

|r„0)|  -  \Cyj  U.l)  (a,  (f)  <7,  (/))  '|  «  1,  (60) 

where  Eq.  (59)  comes  from  Eq.  (B-24),  and  Eq.  (60)  is  obtained  by  analogy  to  Eq.  (59). 

It  follows  from  Eqs.  (47)  and  (49)  that  the  threshold  crossing  rate  for  stochastic  threshold 
functions  Y„(i)  may  be  written  as 

oo 

mj(i)  -  J  fd(di)  Y  ij).  (61) 

—  oo 

Equations  (59)  and  (60)  permit  considerable  simplification  of  the  functions  appearing  inside  the 
integral  in  Eq.  (6i).; 
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From  Eqs.  (55)  and  (59): 


(Ty 


Mu)[4>(p)  +  />4>(p)j, 


(62) 


where  u,  and  <!>(•),  are  given  by  Eqs.  (53),  (57),  and  (58),  respectively.  From  Eqs.  (56) 
and  (59), 


P  “  (m,-i\)laz. 


By  analogy  with  Eq.  (52),  and  making  use  of  Eq.  (60), 


/>  r  - 


•ry.!  4>(u„)  \  hr,'  <#>(p„) l. 


where  u„  and  p„  are  defined  similarly  to  Eqs.  (S3)  and  (63)  as 

u„  s  ({  -  mY) /ar 
and 

P„  »  <«/  -  0/*z  • 


(63) 

(64) 

(65) 

(66) 


It  follows  from  Eqs.  (61)  •  (66)  and  some  algebra  that 


mjd)  -  (2ir)  ,/J 


9vv. 


where,  by  definition 


and 


my  -my 

F(m,-m,). 

(67) 

a  yy 
'  *11 

MW 

n 

(68) 

l  4. 

>  r  ,r  1  / 

(69) 

The  quantities  »»,  (/).  <r »  (/),  and  <r  ,  (/),  in  Eqs.  (67)  -  (69)  are  calculated  by  means  of 

Eqs.  (39),  (40),  and  (41),  respectively.  Finally,  the  expected  number  of  threshold  crossings  in 
a  time  interval  T„  may  be  obtained  as 

iiijiO.T)  -  fiiijU)  Ji,  (70) 


with  »»;(/)  given  by  Eq.  (67). 

Equation  (67)  provides  the  basis  for  analyzing  a  much  broader  range  of  possible  adaptive 
threshold  schemes  than  Fig.  2  might  suggest.  For  example,  straightforward  generalizations  of 
Eq.  (67)  may  be  applied  to  the  structures  depicted  in  Figs.  9-11. 


19 


20 


'*'?HV  <”=*'»**»  **,  *n*i  ^■•'•a-wsw^wrav.-  * 


NRL  REPORT  8320 


Ft*  II  —  The  tapped  delay  line  udaplivc  threshold  scheme 
shown  here  generalizes  the  structure  of  Fig.  10.  The  ir.insfcr 
functions  //„</ ),  »  - 1,2.  .  2\  at  each  of  the  2  V  taps  of  the 

tapped  delay  line  are  design  variables,  chosen  to  maximize  the 
receiver's  performance  ugamst  a  particular  background  scene,  or 
set  of  background  scenes  There  is  no  u  /uion  reason  why  the 
various  delays  should  be  chosen  as  equal  to  one  another,  other 
than  fo'  fabricalionul  simplicity.  More  generally,  additional 
degrees  of  freedom  are  incorporated  by  allowing  these  delays  to 
lake  on  distinct  values. 


While  the  function  F(ni/  -in/)  in  Eq.  (67)  is  fairly  complicated.,  its  zero-ordinate  is 

unity: 

f( 0)  -  1,  (71) 

Equation  (71)  is  a  highly  desireable  feature,  as  explained  in  the  next  section. 


It  is  noted  that  for  uniform  backgrounds,* 


(w,  /tT))1  -  (A,/;,/ A/) 

and 

(<r  ,  /<7  ,)3  -  (A/,/A/)\ 


(72) 

(73) 


where  A/ is  the  noise  bandwidth  of  //(/),  and  A,/„  is  the  noise  bandwidth  of  //„(/)  (cf.  Fig. 
2).  It  follows  from  Eqs.  (68),  (69),  (72),  and  (73)  that 

trj,  =  |l+  (A/;,/A/)J  (74) 

and 


ffi  1+  (A/, ,/A/)3 


(75) 


Eq  (73)  is  domed  by  assuming  a  rectangular-shaped  //(/)  having  an  upper  cut-off  frequency  /„  and  a  noise 
bandwidth  \i  --  /„ 
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Equations  (74)  and  (75)  are  good  approximations  so  long  as  my(i)  is  slowly  varying  compared 
to  //(/)  and  /»„(r), 

DISCUSSION 

Experience  with  the  numerical  evaluation  of  Eq,  (70),  with  >hj(t)  given  by  Eq,  (67),  has 
shown  that  the  principal  contributions  to  the  integral  mj  of  mj  accrue  in  the  neighborhood  of  a 
discrete  set  of  times. 

Moreover,  it  has  been  shown  that  these  important  discrete  times  are  of  two  types: 
"mean-crossing  times,"  and  "times  of  closest  approach"  [)7]. 

Mean-crossing  times  satisfy  the  following  two  conditions  simultaneously:. 


W  y  ( tf)n  )  *  M  y  ^  (  f,tn  ) 

( 76 ) 

(77) 

For  each  solution  of  Eq.  (76)  that  satisfies  constraint  (77),  i.e.,  each  time  the  mean  current 
w,(r)  crosses  the  mean  threshold  w,  (r)  with  positive  slope,  the  expected  number  of  crossings 
nij  is  incremented  by  unity.  Whenever  mean-crossings  exist  during  the  interval  T„f  it  is  gen¬ 
erally  not  necessary  to  perform  the  integral  of  Eq.  (70);  in  this  case,  the  expected  number  of 
threshold  crossings  in,{0,T)  is  well-approximated  by  the  number  of  mean-crossing  limes  /„„ 
during  the  interval  T„.  Clearly,  it  is  desireable  that  no  mean-crossing  times  exist  except  when 
there  is  a  target  in  the  scene,  this  may  be  taken  as  a  reasonable  first  principle  of  search  set 
design  for  operation  against  structured  backgrounds. 

If  »i*(r)  lies  below  »»,  (/)  on  the  lime  interval  i.e.,  if  there  are  no  mean-crossings 
during  T„ ,  the  crossing  count  integral  Eq,  (70)  is  generally  uominated  by  contributions  accruing 
in  the  neighborhood  of  "closest  approach  times"  where  by  definition 

»»/(/,„)  =*  n i/  </,,,)  (78) 

It  follows  from  Eqs.  (67),  (71),  and  (78),  that  the  mean  threshold  crossing  rate  for  adaptive- 
threshold  (AT)  processors  is 

,  ,  •»  >„-'»> 

/«,/(/)  =  (27r)  /2 - - -  (79) 

in  the  neighborhood  of  all  closest  approach  times.  Thus,  the  complicated  function  /■(•)  appear¬ 
ing  in  Eq.  (67)  is  generally  not  needed  in  evaluating  Eq.  (70). 

Analogous  to  Eq.  (79),  the  mean  threshold  crossing  rate  for  fixed-threshold  (FT)  proces¬ 
sors  is  given  by 

,,,  w /  V'„  hi*  I 

/«(,(/)  =  riijii  |  y„,0)  =  (2n)  2  - —  <t> - 1-  (80) 

ir  t  ,r  >  I 
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Equation  (80)  is  derived  from  the  background-conditional  crossing  rate,  Eq.  (62),  by  imposing 
a  saddle  point  condition  analagous  to  Eq.  (78):, 

p  -  (m/  -  yt)lv,  -  0.  (81) 


Alternatively,  Eq.  (80)  may  be  derived  by  noting  that  for  fixed  threshold  processors: 

(tI  —  cr  ?,  =0 

»  'o 


and 


'«v,  =  y„. 


From  Eqs.  (68),  (69),  and  (82), 
and 


<rlY  5  <t] 
«}/„  =  >*/ 


(82) 

(83) 

(84) 

(85) 


for  fixed-threshold  processors.  Substituting  Eqs.  (83)  •  (85)  into  Eq.  (79),  the  adaptive- 
threshold  crossing  rate  >iitr  is  seen  to  reduce  properly  to  the  fixed-threshold  crossing  rate  mir 
given  by  Eq.  (80). 

The  advantages  and  disadvantages  of  adaptive  threshold  (AT)  processing  vis-a-vis  fixed 
threshold  (FT)  processing  may  be  evaluated  by  comparing  the  crossing-rate  expressions  in  Eqs. 
(79)  and  (80), 

It  follows  from  Eq.  (80)  that  fixed  threshold  processors  will  suffer  background-induced 
mean-crossings  whenever  the  peak  target  amplitude  is  less  than  the  clutter  amplitude.  This 
situation  is  depicted  in  Fig.  6.  Thus,  it  may  be  said  that  target-to-dutter  ratios  less  then  unity 
cause  the  FT  processor  performance  to  be  “background-structure-limited"  (BSL),  In  this  case, 
each  "false  alarm"  can  be  associated  with  a  structural  feature  in  the  backgound.  The  effect  of 
quantum  noise  (as  reflected  in  the  magnitude  of  <r )(  for  example)  is  then  totally  overshadowed 
by  background  structure  effects. 

Inspection  of  Eq.  (79)  and  Fig.  7  shows  that  an  Adaptive  Threshold  (AT)  processor  need 
not  suffer  background-induced  mean-crossings.  If  the  filter  //„(/)  (cf.  Fig.  2  )  can  be  chosen 
such  that  /wt  (/)  "tracks"  the  background-induced  variations  in  /«>(/),  the  background-induced 
mean-crossings  can  be  eliminated. 

Such  false  alarms  as  then  occur  are  distributed  randomly  in  time,  and  are  not  associated 
with  particular  features  in  the  background  scene:  the  residual  false  alarms  are  due  to  quantum 
noise.  An  1RST  sensor  operating  in  this  regime  (e  g..  Fig.  7)  is  said  to  be  "quantum-noise- 
limited"  (QNL)  in  its  performance.  Clearly,  QNL  (quantum-noise-limited)  operation  is 
preferable  to  BSL  (background-structure-limited)  operation. 

It  should  be  noted  that  the  adaptive  threshold  performance  advantage  just  described  is 
only  realized  when  the  background  scene  is  non-uniform.  The  performance  of  AT  processors  is 
generally  inferior  to  the  performance  of  FT  processors  when  the  background  scene  is  uniform 
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and  of  known  brightness,  in  this  case,  the  adaptive  threshold  false  alarm  rate  ( FARat )  is 
greater  (i.e.,  worse)  than  the  fixed  threshold  false-aiarm  rate  (FARfT)\ 

FAR^t  ^  FAR/. (86) 

In  order  to  make  the  fals*  alarm  rate  comparison  above  meaningful,  it  is  assumed  that  the 
adaptive-threshold  gain  K  in  Fig.  2  has  been  adjusted  to  achieve  equal  target-detection  sensitivi¬ 
ties  for  the  two  processors  being  compared.* 

A  false-aiarm  penalty  (FAP)  is  now  defined: 

FAP  s  lOlogio  (FARat/FAR/t).,  (87) 

The  false-aiarm  penalty  is  a  measure  of  the  performance  disadvantage  that  accrues  when  an  AT 
processor  is  used  when  it  truly  isn't  needed,  Eq.  (87)  may  be  written  as 

FAP  -  lOlogio  \mATlm/r\  dB.  (88) 

with  mAf  and  given  by  Eqs.  (79)  and  (80),  respectively.  The  evaluation  of  Eq.  (88)  is 
simplified  considerably  by  making  use  of  Eqs.  (74)  and  05).  it  follows  from  Eqs.  (74),  (75), 
(79),  (80),  and  (88)  that: 

my -mr  |: 

FAP  -  2.17  — - +  Slog io  ll-a(l-a)].  dB.  (89) 

I  +a  or  y  | 

where,  by  definition. 


Equation  (89)  simplifies  still  further  in  the  limit  u«l.  Thus 


FAP  -  2.17  « 


my  -my 


dB  ,  a  « 1 , 


As  an  example,  it  is  assumed  that  the  threshold-constant  K  in  Fig.  2  is  adjusted  until  "  the 
threshold  is  five  sigmas  above  the  mean,"  i.e., 

|  my -my  I 


when  the  search  set  is  observing  a  uniform  scene  of  known  brightness.  Furthermore,  it  is 
assumed  that  the  noise  bandwidth  of  the  target  filter  //(/)  in  Fig.  2  is  twice  as  large  as  the 
noise  bandwidth  of  the  threshold-setling-filter  //„(/), 


a-\TfrT 


ll  is  assumed  lii.it  the  threshold  litter  //„(/)  responds  too  slowly  to  suppress  last  rise-lime  target-induced  threshold 
crossings 
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Jt  follows  from  Eqs.  (89),  (92),  and  (93)  that  the  false  alarm  penalty  is  : 

FAP  -  17.46  dB,  (94) 

corresponding  to  a  value  of  FARat  (  adaptive  threshold  false-alarm  rate)  about  56  times  worse 
than  FARft  (fixed-threshold  false-alarm  rate). 

Generally,  A/„  should  be  chosen  smaller  than  A/  to  both  minimize  the  false  alarm 
penalty,  (Eq.  (89)),  and  to  prevent  a  too-rapid  threshold  response  that  would  tend  to  suppress 
target-induced  threshold  crossings.  On  the  other  hand,  A/„  should  be  chosen  large  enough  to 
allow  the  threshold  to  accurately  follow  most  of  the  structure  in  the  background  scene.  Clearly, 
the  choice  for  A/„  involves  degrading  system  performance  against  uniform  backgrounds  for  the 
sake  of  improved  performance  against  non-uniform  backgrounds. 

It  appears  likely  that  a  more  favorable  trade-off  could  be  achieved  with  the  receiver  struc¬ 
ture  shown  in  Fig.  9,  both  from  the  standpoint  of  (a)  decreasing  the  false  alarm  penalty,  Eq. 
(89),  and  (b)  improving  the  background  tracking  properties  of  0).  Equation  (79)  and  the 
entire  analysis  of  the  preceding  section  is  easily  adapted  to  the  structures  of  Figs.  9-11.  The 
false-alarm  penalty,  Eq.  (89),  decreases  roughly  as  (2\)',/2  for  the  detector  of  Fig.  10.  The 
improvement  in  uniform  background  performance  thus  obtained  for  large  values  of  N  is  gained 
at  the  expense  of  degraded  performance  against  cluttered  scenes,  as  compared  to  detectors  with 
small  values  of  N,  The  good  background-tracking  capability  of  the  structure  in  Fig.  9  combined 
with  the  low  false-alarm  penalty  of  the  structure  of  Fig.  10  can  be  obtained  by  employing  a 
two-dimensional-detector  array  with  time-delay  and  integration  (TDD  logic. 

In  order  to  put  this  discussion  on  a  concrete  quantitative  basis,  particular  background  and 
target  radiance  distributions  must  be  chosen,  and  the  mean  current  /«»(/)  calculated  by  means 
of  Eq.  (A-2)  in  Appendix  A.  The  target  detection  and  clutter  rejection  capabilities  of  a  given 
candidate  adaptive-threshold  processor  can  then  be  analyzed  by  means  of  Eq.  (79),  Intercom¬ 
parisons  of  the  numerical  results  thus  obtained  for  a  variety  of  different  processor  structures 
should  then  allow  quantitative  conclusions  to  be  drawn  concerning  such  issues  as: 

•  The  performance  penalty  caused  by  failing  to  match  the  sensor's  instantaneous  field- 
of-view  to  the  angular  size  of  the  target. 

•  The  potential  performance  advantages  of  time-delay  and  integration  (TDD, 

•  The  best  value  of  N ,  and  the  desireability  of  having  different  transfer  functions  //„(/) 
for  each  of  the  2  A  taps  in  the  tapped  delay  line  structure  of  Fig.  11. 

•  The  advantages  that  may  be  gained  by  employing  two-dimensional  threshold  process¬ 
ing,  in  which  the  "target  signal"  Y(i)  and  threshold  function  >',,(/)  are  derived  from  detectors 
scanning  at  different  elevations. 

The  only  important  obstacle  to  performing  analyses  of  the  kind  described  above  is  the  lack 
of  high-spatial-resolution,  radiometric,  infrared  background  imagery. 
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Appendix  A 

CALCULATING  THE  AVERAGE  PHOTOCURRENT 
FROM  BACKGROUND  DATA 

The  objective  of  this  Appendix  is  to  present  Eq.  (A-2),  which  expresses  the  average  value 
of  the  random  current  XU)  (cf.  Fig.  1), 

mxU)  -  £U0)},  (A-l) 

as  a  function  of  the  scene  radiance  distribution  mLU)  and  a  number  of  important  sensor  paranv 
eters. 

As  derived  in  Ref.  17, 

mxU)  -  K„R,  J  MTF(f)  9(0  my(t)  exp(/2irf  •  v/)</f.  (A-2) 

The  various  quantities  appearing  on  the  right-hand-side  of  Eq.  (A-2)  will  now  be  defined. 

The  constant  K„  is  given  by 

Ku-*rJ(2f*)1.  (A-3) 

where  r„  is  the  transmittance  cf  the  optics,  and  ./*  is  the  focal  length  ratio  of  the  optics. 

The  constant  Rt  is  the  current  responsivity,  given  by: 

Ri  m  *if/h»,  (A-4) 

where  y  is  the  detector  quantum  efficiency,  e  is  the  electronic  charge,  li  is  Planck's  constant, 
and  v  is  the  average  optical  frequency  of  the  incident  light.: 

The  variable  of  integration  in  Eq.  (A-2),  f,  is  the  two-dimensional  vector  spatial  fre¬ 
quency.  The  quantity  MTF(f)  is  the  modulation  transfer  function  that  characterizes  the  image 
blurring  effect  of  the  optical  train,  normalized  such  that 

MTF(O)  -  1.  (A-5) 


For  the  detector  geometry  depicted  in  Fig.  Al,  the  quantity  Pit)  is  defined  as  follows: 


/Mr) 


1  r  €  Cf,*,, 

0 


(A-6) 
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Fig  Al  -  Focal  Plane  Geometry.  The  focal 
plane  irradiance  is  stationary  in  the  coordi¬ 
nate  system  with  origin  ft  Vector  r„  locales 
the  center  of  a  detector  of  area  /t^,,.  For 
scanning  sensors.  r0  is  a  function  of  lime 


The  function  fit )  is  defined  as  the  two-dimensional  spatial  Fourier  transform  of  /Mr): 

fit)  -  J  Pit)  exp(./2irfr)</r,  (A-7) 

Similarly,  wy(f)  is  the  two-dimensional  spatial  Fourier  transform  of  the  radiance  distribution 
m,  it),  It  should  be  noted  that  in,  it),  like  «i»(/)  in  Eq.  (A-2),  is  an  ensemble  average  value 
over  the  photon  fluctuation  statistics  of  the  incident  light. 

The  radiance  distribution  in,  (r)  is  characteristic  of  a  particular  infrared  scene,  and  may  be 
estimated  by  means  of  a  radiometric  Thermal  Imaging  System  (1,21  (TIS)  of  higher  spatial 
resolution  than  the  model  system.  It  is  also  highly  desirable  that  the  dwell  lime  of  the  TIS  be 
much  longer  than  that  of  the  model  system,  since  the  analysis  requires  knowledge  of  the  mean 
radiance  of  the  scene  established  by  averaging  over  the  photon  fluctuation  statistics  of  the 
incident  light*. 

The  spectral  filter  chosen  for  use  with  the  TIS  should  match  the  combined  TIS  optical 
train/photodetector  spectral  response  to  that  of  the  model  system.  This  is  necessary  because 
there  is  no  way  to  reliably  calculate  the  radiance  of  a  scene  measured  in  a  waveband  in 
terms  of  the  radiance  of  the  same  scene  measured  in  a  different  spectral  band  AA;. 

The  quantity  v  in  Eq.  (A-2)  is  the  focal  plane  scan  velocity,  which  may  be  calculated  as 

|v|  =  2W*  DJT,.  (A-8) 

where  DUI,  is  the  diameter  of  the  optical  aperture,  and  T,  is  the  system  frame  lime. 


'Averaging  out  the  photon  noise  by  either  increasing  the  TIS  dwell  lime  or  b>  performing  Ir.une  .iilihlion  decreases  the 
noisiness  ol  the  fIS  image,  and  improves  the  goodness  ol  the  TIS  imager',  as  an  estimate  ol  m,  <r> 
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Appendix  B 

NOISE  CURRENT  CORRELATION  FUNCTIONS 

The  object  of  this  Appendix  is  to  derive  Eqs.  (BOO)  •  (B04),  which  are  needed  to  evalu* 
ate  Eq.  (55)  for  the  threshold  crossing  rate  thjU\y„.  yj. 

The  starting  point  for  this  discussion  is  the  linear  system  input/output  relation  between 
the  random  processes  X(i)  and  Y(i)  (cf.  Figure  1): 


Y(i)  -  f  </m  X (/i) Ii(i— it), 

ao 

It  follows  directly  from  Eq.  (B-l)  that 

d(i<lK  C\\(k,n)  ll(l[—K) 

where  the  covariances  <’>>  and  (\\  are  defined  by 

<•,,</„  /2)  s  A’(( K(/,)  -/«»</,))  |  K</2)-/M> </j)|)„ 
<  h(A.m)  s  /;||TU)-w,(A)I 

and  where 

m,0)  s  Al  K</)}.- 


and 

With  the  definition 
it  follows  from  Eq.  (B-3)  that 

and 


«»»(/)  =  fc’U(r)}. 

7Ai)  =  Yii), 

(  fj)  =  d,,  |(  > > ( r | ,  / j) ) 

(■'//ill-  h)  “  9/,  It  i»(*i>  '2>l- 


From  Eqs.  (B-2),  (B-8),  and  (B-9), 

gt  gt  no 

C)/(t l})  =*  J J  ilyulk  C\\(K,,n)  Ii(i |  -  A)  d,2h(t2-n) 
and 

v  '//(>[.  I2)  =  J* J  a  (Jfulk  C\\(K,n)  -  A)  d,,/i(/ 2~m)- 

Setting  / 1  =  / 1 = /  in  Eqs.  (B- 10)  and  ( B- 1 1 ) ,  and  noting  that 

C//U.,  l)  =  «r/(/)  =  E[Z2(l)  |  -  «//(/),; 


(B-l) 

(B-2) 

(B-3) 

(B-4> 

(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 

( B- 1 0) 

(B-l  I) 

(B-12) 


28 


NRL  REPORT  8320 


it  follows  that 


and 


Cy/(i,  I )  ■  JJ  ^  dfidk  Cxx(k,  (i)  h (t  ~ k)  &,h(t -k) 

cr|(f)  -  //"  dfidk  Cxx(k.fi)  d ,hU-k)  B.h(l-fi). 

Setting  /)  -  f2“  i  *n  Eq.  (B-2),  and  noting  that 

CYy(i,i)  -  <r\(i)  -  E\rHt)\  - 


it  follows  that 


<t\ U)  -  f  f_m  dfidk  Cxx(k.  fi)h(i  -  k)h(t-fi).: 


(B-li) 

(B-14) 

(B-15) 

(B-16) 


An  expression  for  Cxx(k,n )  is  now  required  before  the  analysis  can  be  carried  any 
further.  An  adaptation  of  Eq.  (4.3.13)  on  p.  115  of  Ref.  28  leads  to: 

cr\(t)  ~  e  f  ^  dn  lt2(t -fi)  (B-17) 

where  e  is  the  electronic  charge.  Consistency  between  Eqs.  (B-16)  and  (B-17)  requires  that 

C \y(k.  fi)  —  e  ntx(fi)  b(k-fi).  (B-18) 

where  6(.)  is  the  Dirac  deln  function.  Covariance  functions  like  Eq.  (B-18)  are  characteristic 
of  non-stationary  white  noise  (26,29).  From  Eqs.  (B-13),  (B-14),  and  (B-18), 

(B- 19) 
(B-20) 

(B-21) 


and 


Noting  that 


CyzU.  t)  »  e  I  dfi  li(i-fi)  B,h(t-fi)  my(fi) 

•r  —  OO 

•  ef_„  dfi  [B,h(i  -jt)P  mx(fi)( 


h(t  -fi)  B,li(i  -fi)  *  —  B,h2(i  -fi), 


it  follows  from  Eqs.  (B-19)  and  (B-21)  that 

Cy/(l.  i)  -  ~  dfi  h2(i-fi)  Itty(fl) 

From  Eqs.  (B-17)  and  (B-22), 

1  a 


Defining  r(t)  as 


(.’>/(/,  I)  -  —  d,!o-|(/’»j  "  <Ty(l)  a  >•(/), 


r(f)  =  C)z(l,  I )  [<T  y(t)<T/(t)] 


(B-221 


(B-23) 


(B-24) 
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it  follows  from  (B-23)  and  (B-24)  that 

r(t)  -  [&  y(t)/(TZ(t)). 

(B-25) 

Taking  the  expected  value  of  both  sides  of  Eq.  (B-l)  leads  to  the  result: 

niy(i)  -  J  m  dfi  my(ti)  hit  -fi). 

(B-26) 

Taking  the  time  derivative  of  (B-l),  we  have 

Yd)  -  ZO)  -  dfl  xift) 

(B-27) 

Taking  the  expected  value  of  both  sides  of  (B-27),  we  have 

/»/(/)  -  J  ^  dn  m.y(n)  B,h(t  ~n).. 

That  is, 

|  \ 

m/U )  -  0,1/.,,,  dfi  »Mm)  It(i-n)  J. 

(B-28) 

From  Eqs.  (B-26)  and  (B-28), 

mz(t)  -  thy(l). 

(B-29) 

Defining  the  convolution  operator  as  in  Eq.  (37a),  Eqs.  (B-26),  (B-17),  (B-20), 
(B-25)  may  be  written  as: 

(B-29),  and 

myO )  -  /»(/)  ®  my(t). 

(B-30) 

a-](i)  -  e  li2(i)  ®  mt0), 

(B-31) 

<r2/(t)  -  e[li(l)]2  ®  /»*(/), 

(B-32) 

m/(i)  m  mt(f). 

(B-33) 

and 

r(l)  -  {a  t(l)/a/(l)}. 

(B-34) 

Equations  (B-30)  •  (B-34)  are  the  desired  results. 
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Appendix  C 

RELATIONSHIPS  BETWEEN  FAR,  PD,  and  ntj 
Relationship  Between  FAR  and  mj 

The  complete  description  of  an  IRST  sensor's  performance  under  a  given  set  of  opera- 
tional  conditions  requires  the  simultaneous  specification  of  both  the  False  Alarm  Rate  (FAR) 
and  the  Probability  of  Detection  (/*/>)  for  a  "target"  within  the  sensor’s  field  of  view.  However, 
both  "false  alarms"  and  target  detections  are  manifested  as  threshold  crossings  by  the  signal  pro¬ 
cessor.  Thus,  the  object  of  this  Appendix  is  to  relate  the  traditional  IRST  performance  meas¬ 
ures,  Pn  and  FAR,  to  the  expected  number  of  threshold  crossings  n?j  over  prescribed  intervals 
of  time, 

it  is  assumed  that  the  average  current  m*(r)  is  known*  on  an  interval  of  time  |/|  <  772, 

The  expected  number  of  threshold  crossings  on  the  interval  |/|<772  is  defined  as 
mj( 0,  T),  Defining  the  false  alarm  rate  as  the  expected  number  of  threshold  crossings  per 
"reference  interval"  T,sf,  the  following  relationship  obtains  between  FAR  and  my 

FAR  -  (Tttf/T)  77,  (C-l) 

For  example,  if  FAR  is  defined  as  the  average  number  of  false  alarms  per  week,  T,t(  is  set 
equal  to  the  number  of  seconds  in  one  week;  if  FAR  is  defined  as  the  average  number  of  false 
alarms  per  system  dwell  lime,  then  7ref  is  set  equal  to  the  dwelt  time  (again  expressed  in  units 
of  seconds). 

It  is  implicitly  assumed  in  Eq.  (C-l)  that  the  scene  under  observation  does  not  include  a 
target,  so  that  each  threshold  crossing  that  occurs  gives  rise  to  a  "false  alarm," 

Relationship  Between  P,)  and  my  First-Order  Approximation 

Although  not  as  straight-forward  as  Eq.  (C-l),  a  relationship  between  Pt)  and  mj  can  also 
be  established. 

As  prelude  to  the  definition  of  P„,  a  "decision  interval"  T»  is  first  defined.  The  interval 
T/t  is  presumed  to  bracket  the  entire  period  of  time  during  which  'he  current  Y(i)  manifests 
target-induced  fluctuations. 

Assuming  that  a  target  is  present  in  the  scene,  the  number  of  threshold  crossings  that 
occur  during  the  interval  Tt)  is  defined  as  the  integer  random  variable  J.  The  discrete  probabil¬ 
ity  density  function  of  J  is  denoted  as 


The  quantity  »iv(/l  may  he  specified  u  ;«nun,  or  it  may  be  calculated  in  terms  of  the  radiance  of  a  particular  back¬ 
ground  scene  (as  discussed  in  Appendix  A) 
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The  probability  of  detection  Pp  is  now  defined  as  the  probability  that  one  or  more  thres¬ 
hold  crossings  occur  during  the  decision  interval:. 

Pp-i  M  (C-2) 

/-l 

(The  likelihood  of  a  background-induced  crossing  during  Tp  has  been  neglected.)  Unfor¬ 
tunately,  the  problem  of  obtaining  a  formulation  for  fj  appears  to  be  quite  difficult  [30] .  The 
focus  of  this  paper  has  been  on  the  development  of  formulations  for  the  expected  number  of 
threshold  crossings: 

ntj  a  £{./}-  £y/,0),  (C-3) 

/-i 

In  order  to  establish  a  relationship  between  Pp  and  mj,  Eqs.  (C-2)  and  (C-3)  are  written  as: 

Pp-fM)  +  i/jW  (C-4) 

(-2 

and 

mj-fj(\)  +  ZJfjii).  (C-5) 

i-l 

Assuming  that  the  probability  of  two  or  more  threshold  crossings  is  negligible  during  the  deci¬ 
sion  interval  Tp,  Eqs.  (C-4)  and  (C-S)  can  be  estimated  as 

Pp  ~M\)  (C-6) 

and 

=*/,(!).  (C-7) 

It  follows  from  (C-6)  and  (C-7)  that 

Pp  =  mj(0.  T,,)..  (C-8) 

According  to  Eq.  (C-8),  the  expected  number  of  threshold  crossings  during  the  decision  inter¬ 
val  Tp  provides  a  good  estimate  of  the  detection  probability  Pn ,  so  long  as  the  probability  of 
two  or  more  crossings  during  Tp  is  negligible. 

Relationship  Between  Pp  and  /«,:  Pr  posal  for  a  Second-Order  Approximation 

It  appears  only  reasonable  that  an  improved  estimate  could  be  obtained  for  Pp  if  the  vari¬ 
ance  <r]  of  J  were  known  in  addition  to  the  mean  mj  of  /  It  will  now  be  shown  how 
knowledge  of  a  j  can  be  used  to  calculate  a  second-order  approximation  for  Pp  (compare  with 
Eqs.  (C-4)  and  (C-6)): 

Pn  =  fj(\)  +  Jj(2).  (C-9) 

Unfortunately,  the  technique  described  in  this  section  for  calculating  the  second-order 
approximation  to  Pp  cannot  be  implemented  until  a  formulation  for  a)  is  developed  analagous 
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to  Eq.  (67)  for  rhj.  In  this  connection,  it  is  noted  that  Bendat  has  derived  an  equation  for  the 
crossing  count  variance  of  stationary  processes  (31)  His  result  (cf.  also  Ref.  30)  is  far  more 
complicated  than  the  analagous  Eqs.  (8)  and  (9)  for  mj.  Thus,  a  generalization  of  Bendat's 
result  for  crj  to  the  case  of  nonstationary  processes  and  stochastic  threshold  functions  may 
prove  to  be  a  difficult  problem.  Nonetheless,  it  is  now  assumed  that  a  formulation  for  crj  can 
be  obtained,  analagous  to  the  development  for  m.,  as  noted  above  and  found  in  the  main  text, 
Analagous  to  Eq.  (C-3), 


crj  = 


/-I 


-*»/• 


(C-10) 


Substituting  Eq.  (C-3)  into  Eq.  (C-10)  leads  to  the  following  expression  for  crj: 

crj  -  Ml)  [l-/;(l)|  +  4/, (2)  [l-/y(l)-/,(2)j  +  Ej.  (C-U) 


where 

Ej  ■  £yJ/-2  1/(1)  +  2/(2)|  ZJfAiA.  (C-12) 

/-}  1  1  ,-J  (;-J  j 

Assuming  that  (he  probability  of  three  or  more  threshold  crossings  is  negligible  during  the  deci¬ 
sion  interval,  Eqs.  (C-3)  and  (C-tl)  are  approximated  as: 

"j  =*  MD  +  2/,(2)  (C-13) 

crj  a*  Ml)  [l  -/,(!)(  +  4//(2)  [l-//(D  -//(2)|.:  (C-14) 


Calculation  of  mt  from  Eqs.  (67)  and  (70),  and  an  analagous  calculation  for  a  j,  enables  Eqs. 
(C-13)  and  (C-14)  to  be  solved  for  approximations  to  ./>(1)  and  /,(2).  The  second-order 
approximation  to  Pl}  is  then  obtained  by  means  of  Eq.  (C-9). 

If  Eq.  (C-9)  is  found  to  yield  an  appreciably  different  result  than  Eq.  (C-6),  third-order  or 
even  higher-order  approximations  to  P„  may  be  required;  otherwise,  the  first-order  Eq.  (C-6) 
is  then  verified  as  a  good  approximation  for  P». 
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